NAS A Technical Memorandum 106642 




t 

k 

Theoretical, Experimental, and Computational 
Evaluation of a Tunnel Ladder 
Slow-Wave Structure 


Thomas M. Wallett and A. Haq Qureshi 
Lewis Research Center 
Cleveland, Ohio 


June 1994 



National Aeronautics and 
Space Administration 


THEORETICAL, EXPERIMENTAL AND COMPUTATIONAL EVALUATION 
OF A TUNNEL LADDER SLOW-WAVE STRUCTURE 


Thomas Michael Wallett and A. Haq Qureshi* 
National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


SUMMARY 

The dispersion characteristics of a tunnel ladder circuit in a ridged wave guide were 
experimentally measured and determined by computer simulation using the electromagnetic code 
MAFIA. To qualitatively estimate interaction impedances, resonance frequency shifts due to a 
perturbing dielectric rod along the axis were also measured indicating the axial electric field strength. 
A theoretical modeling of the electric and magnetic fields in the tunnel area was also done. 


INTRODUCTION 

Of a variety of structures considered in recent years for possible applications at millimeter- 
wave frequencies, one that has received attention is the thin ladder circuit which was presented 
conceptually by Karp in 1955 [1] and its modification discussed by him in 1960 [2]. Similar parallel 
line and ladder structures were analyzed by Pierce [3] and Butcher [4] without space harmonics and by 
Froom et al. [5] and Kosmahl and O’Malley [6] using space harmonics. The predicted dispersion 
characteristics were in close agreement in all, however the values of interaction impedance varied. In 
their presentation, Kosmahl and O’Malley suggested that high beam operating voltages in the range of 
tens of kilovolts would make possible a non-space harmonic operation of a forward-wave ladder based 
amplifier at millimeter wavelengths. They also suggested that the gain rate would further benefit from 
the high interaction impedance associated with the ladder. Kosmahl and Palmer conducted an analysis 
of the idealized or modified Karp circuit referred to as the TunneLadder [7]. Their results pointed out 
to the suitability of this structure as a high impedance, narrow bandwidth circuit of about 1% that is 
voltage tunable over a frequency range of 5% and has excellent heat handling capability. 



Figure 1 - Idealized circular tunnel ladder structure. 


*National Research Council — NASA Research Associate at Lewis Research Center. 


A simplified theoretical model of the TunneLadder structure was used by Kosmahl to predict 
the dispersion characteristics and interaction impedances (Figure 1). This model uses a circular wire 
loop to represent the hexagonal cross section of the actual tunnel region. The electric and magnetic 
fields are represented in cylindrical geometry inside the loop region and by Cartesian geometry outside 
• the region. The boundary conditions are then applied at the various regions. 

The actual TunneLadder is a novel, high impedance structure that was introduced by Karp 
and does not lend itself to a closed form analysis. The idealized structure, although quite simplified, is 
more amenable to closed form solution. With the idealized tunnel cross-sectional area equal to the 
actual cross-sectional area, the model predicts at least qualitatively the dispersion characteristics of 
both the symmetric and antisymmetric ladder modes. 

A 29 GHz TWT utilizing the TunneLadder structure was fabricated and tested by Varian 
under a NASA contract [8,9]. The test results confirmed the potential applicability of the ladder-based 
structures for use in millimeter-wave traveling-wave tubes. 


THEORETICAL EVALUATION OF THE ELLIPTICAL TUNNEL 

A realistic approach would be to approximate the actual hexagonal tunnel cross section as an 
ellipse supported by dielectric slabs (Figure 2). 



Figure 2 - Idealized elliptical tunnel ladder structure supported by dielectric slabs. 

For this analytical effort, elliptic cylinder coordinates would be the preferred choice for field 
% representation. The result would be a model to closely represent the electric and magnetic fields inside 

the tunnel region. 

• Electromagnetic waves propagating in elliptical metal pipes have been investigated as early as 

1938 [10]. Since then several papers dealing with wave propagation [11], computation of critical 
frequencies [12] and cutoff wavelengths [13,14,15], and attenuation characteristics for elliptical [16,17] 
and dielectric-tube wave guides [18] have been published. 
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Cartesian coordinates in three dimensions are related to elliptic cylindrical coordinates by 
x = d cos(^) cosh(/i), y = d sin(i?) sinh(/i), and z — z where d is half the interfocal distance. The z-axis 
is chosen to be along the direction of propagation with the axial propagation constant designated as /?. 
Propagating solutions in the z-direction are of the form The remaining transverse Helmholtz 

equation V 2 ^ + k 2 xp = 0 where k 2 = (w/c) 2 - /? 2 in the two-dimensional elliptic coordinate system is 

+ d 2 k 2 T [cosh\ti) “ cos 2 (tf)] = 0. (1) 

Ofi ov 

Assuming xp = Af(/i) 0(t?), the transverse Helmholtz equation separates into two equations with 
coupled eigenvalues, h = k T d and b: 

d 2 M/dfi 2 + [ h 2 cosh 2 (/i) — b]M = 0 (2) 

and 

d 2 0/d 2 tf 4 [b-h 2 cos 2 (t?)]0 = 0. (3) 

Solutions for 0 are the Mathieu functions Se m and So m while the corresponding radial 
solutions for M are the modified Mathieu functions of the first kind, Je m and Jo m , and of the second 
kind, Ne m and No m , where m refers to the particular order of the function [19]. The functions Se m (h, 
cos are even in i? with eigenvalues b designated as be m (h) while the functions So m (h, cos tf) are odd 
in i9 with eigenvalues b designated as bo m (h). For a given h, the eigenvalues are ordered such that 
be 0 < bo 1 < be t < ... < bo m < be m < bo m + 1 .... 

Inside the elliptical tunnel the continuity of xp in value and slope at \i = 0 prohibits 
corresponding radial solutions of the second kind, 7Ve m and ATo m . Therefore when the solution for © is 
5e m (/i, cos tf) then the corresponding radial solution for M becomes Je m (h, cosh fi). Likewise, when 
the solution for 0 is So m (/i, cos tf) then the corresponding radial solution for M becomes Jo m (h, cosh 
/i). Also on the metallic ellipse boundary at [i = /i 0 , the z-component of the electric field must be zero. 

There exist transverse magnetic (TM) and transverse electric (TE) modes in an elliptical wave 
guide not to mention hybrid modes. The modes are designated TMj^ or TEj^ where the index / is 
either c (standing for cosine or even) or s (standing for sine or odd). The m is an integer representative 
of the order, and n represents the nth parametric zero of the Mathieu function. For tube devices only 
the TM modes are useful. The field configuration for the lowest transverse magnetic (TM) mode in an 
elliptical wave guide is the TM c01 mode; the TM 501 does not exist [20]. The fields for this mode are 
well known and given in terms of the zero-order, even Mathieu functions and their derivatives. 


E lt = ( ■ ,jf —) c oi Je 'o ( h ’ cosh P) Se o( h > cos exp{j(u>< - P cQX z)} (4) 

\ lk c c01 / 

E# = ( 1 Coi Je o( h > cosh M) Se' 0 (h, cos tf) exp {j(ut - P c0X z)} (5) 

V lk c cOi J 

E z = C 01 Je 0 (h, cosh p) Se 0 (h, cos tf) exp {j(ut - fi c01 z)} (6) 

E /I = f J?* 0 -) Coi Je o( h > cosh M) Se' Q (h, cos tf) exp { j(w< - P cQX z)} (7) 

y k c coi J 

H d = ( ) C 01 Je' 0 {h, cosh p) Se 0 (h , cos t?) exp {j(ut - P c01 z)} (8) 

\‘*c cOl J 

H z = 0 (9) 
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Here / = 




cosh(2/i) — cos(2tf) 


propagation constant is /3 C c01 


k c c01 is the transverse cutoff wave number, the cutoff axial 
= y\(u 2 /c 2 ) — k 2 c01 , and C Q1 is a normalization coefficient. 


The power flow through the tunnel region can be found from 

Vo 27 r 

P = \\ | (. E ^ - E 0 H*Jl 2 M d/i. (10) 

0 0 

To completely specify the fields for the TM c01 mode, the transverse cutoff wave number k c c01 
must be found. Since this is related to the eigenvalue /i, we must determine h such that E z is zero at 
the ellipse boundary /i = /i 0 . 

Since the eigenvalues h and b depend on each other, they must be found simultaneously. 
Solution of the eigenvalues h and b involves the evaluation of continued fractions [21]. For m an 
integer, the eigenvalues are related through 

(11) 

where the a coefficients are to be determined. The ratios aj/a 0 and a_ 1 /a 0 are given by the continued 
fractions 

a » — ~ ( 12 ) 

a " - 1 16(m/2 + n) 2 + 2 h 2 - 46 + h 2 (a n + i/a n ) 

and 

a » — ~^ 2 . (13) 

a " + 1 16(m/2 + n) 2 + 2 h 2 - 46 + h 2 (a n _ J a n ) 


For the lowest order mode, m = 0, the ratio Oj/oq = a _ 1 /a Q . 

For large positive », f ” ^ ”■ V “ 16(„ 7 2 + ^ 

Substituting these limiting values in equations (12) and (13), the ratios a^/a Q and a_j/a 0 are 
found and the eigenvalue relation in continued fraction form becomes 


26 = h 2 — - 


16(1) 2 + 2/i 2 — 46 — ■ 


16(2) 2 + 2h 2 - 46 - 


16(3) 2 + 2/r — 46 — ... 


(14) 


The condition that E z = 0 at \i = n Q is now needed to finally solve for both h and 6. 
to finding h which satisfies 

Je Q (h, cosh n 0 ) = 0. 


This is equivalent 


(15) 


The Mathieu function Je 0 (h, cosh /x 0 ) can be calculated as a summation of Bessel functions 

Je Q (h, cosh /i 0 ) = 4^2 £ ( - l) n B 2 n J 2n( / ‘ cosh ^o) 

n = 0 


where 



(16) 

(17) 
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The procedure then to solve for the proper eigenvalues, given the particular tunnel dimensions, 
is to pick a value for h and solve for b using equation (14). Using h and b one can find the coefficients 
a n to determine the B 2n coefficients. After solving equation (15) for h using equation (16), the 
procedure is repeated. Once the eigenvalues h and 6 are found then values for k c c01 and /? c c01 are 
calculated. With this the unnormalized fields for the elliptical tunnel region are determined from 
equation (4) through (9). 

A BASIC program has been written to calculate the eigenvalue b for a particular eigenvalue h. 
The actual solution of both eigenvalues simultaneously given a set of tunnel dimensions was achieved 
using the program MATHCAD. To determine the dispersion characteristics of the tunnel ladder 
structure the fields exterior to the tunnel region must be found and only by matching admittances at 
the elliptical boundary can a closed form dispersion relation be realized. Since the region exterior lends 
itself to Cartesian coordinates, only a crude match can be achieved. A relatively simple closed form 
solution of such a complicated problem is not warranted since results can only be qualitative in nature. 
However, if the power and dispersion characteristics of the particular tunnel ladder can be found by 
other means, the use of the E z field for calculation of interaction impedance may prove useful. 


EXPERIMENTAL EVALUATION 

For the experimental determination of the dispersion characteristics a 29 GHz tunnel ladder 
structure similar to the Varian structure has been fabricated (Figure 3). 



Figure 3 - Cross-section of the experimental tunnel ladder structure. 


In fact, two copper ladders manufactured by Varian for contract number NAS3-23347 each 
having a slight depression in the middle of the rungs were placed together to form the tunnel ladder. 
The tunnel was hexagonal in shape with a width of 0.0386 inch and a height of 0.0252 inch. Ladder 
and structure dimensions were measured using an optical microscope. Each ladder consisted of 101 
rungs and was 0.0025 inch thick with 0.0066 inch wide rungs and 0.006 inch wide by 0.1105 inch long 
slot spaces. 

A ridged wave guide was fabricated using electrical discharge machining (EDM) from copper 
block in two symmetrical halves. The width and height of the entire cavity of the ridged wave guide 
were 0.212 inch and 0.098 inch, respectively. The cross-section dimension of the ridges was 
approximately 0.0575 inch wide by 0.0214 inch high. Two amorphous boron nitride (c r = 3.973) 
dielectric slabs, 0.0575 inch wide by 0.0125 inch high and glued to both the upper and lower ridges, 
were used to support both ladders. The ladders were positioned between the two halves of the ridged 
wave guide such that the rungs on both ladders were aligned. Screws were used to finally secure both 
halves of the ridged wave guide. A 17-period length of the test structure was cut using wire EDM. 
Minor dimensional inconsistencies in the structure are apparent and are due to the specific machining 
processes and fabrication methods. 
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Since the anticipated high frequency cutoff was near 30 GHz, the tunnel ladder structure was 
initially tested using an HP 8510C network analyzer in WR-28 wave guide. The signal was coupled 
through small holes drilled into the wide faces of the input and output wave guides at the ends of the 
tunnel ladder structure. The attenuation of the signal through the coupling holes was too great, 

• however, and no resonances were observed. Any attempt to enlarge the holes would diminish the 
shorting effect at the ends of the structure. 

• It was decided to use another HP 8510C network analyzer which can operate in coaxial cable 
up to 40 GHz. A test fixture was machined out of brass in which K-connectors were installed in both 
halves such that the signal coupling was achieved from coaxial probes through small holes drilled in 
each half of the fixture. Grooves to fit the length and width of the entire structure were milled on the 
insides of each half of the fixture to center the tunnel horizontally and vertically with the probes. An 
alignment pin connected both halves (Figure 4). 



Figure 4 - Experimental tunnel ladder structure and test fixture. 


The assembled 17-period tunnel ladder structure and test fixture were connected to the network 
analyzer using 2.4 mm coaxial cable (Figure 5). 


• Figure 5 - Assembled experimental tunnel ladder structure and test fixture. 

Measurements from 10 GHz to 35 GHz of the resonance frequencies for the structure were 
taken with (Figure 6) and without (Figure 7) an alumina rod ground to a 0.025 inch diameter and 
inserted in the tunnel region. 
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Six resonances were observed for both sets of measurements. The K-connectors will exhibit a 
resonance just above 30 GHz, but this offered no problem since transmission through the structure was 
negligible near this frequency. 




lin MQG 
-100.0 mUnit.* 

50. 0 



START 1 0.000000000 GHz 
STOP 35 .000000000 GHz 

Figure 6 - Transmission measurement of the tunnel ladder structure without the dielectric rod. 



lin MQG 
-200.0 nUhits 
100.0 mUnit*/ 



START 10 . 000000000 GHz 
STOP 35 .000000000 GHz 

Figure 7 - Transmission measurement of the tunnel ladder structure with the dielectric rod. 
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COMPUTATIONAL EVALUATION 


The tunnel ladder structure was also modeled using the MAFIA code. The MAFIA code 
allows the generation of elliptical cylinders very easily. One period of the actual hexagonal tunnel 
ladder was therefore modeled as an elliptical tunnel ladder, having the same cross-sectional area with 
the same ladder thickness. The height and width of the hexagonal tunnel were used to find the 
eccentricity e of the ellipse with semimajor axis r a and semiminor axis r b . 


y (width ) 2 — (height ) 2 
height 



(18) 


The semimajor and semiminor axes were then found simultaneously from the tunnel area using 

nr a r b = Area ( 19 ) 

The dimensions of the simulated ridged wave guide and dielectric supports were chosen as close 
to the actual dimensions as possible such that bilateral symmetry was preserved vertically and 
horizontally with respect to the center of the tunnel. 

A simulated dielectric rod was given the same dimensions as the actual alumina rod and was 
located in the elliptical tunnel. The shape of the simulated rod cross-section is not quite circular. This 
is due to the resolution of the mesh. A finer resolution would improve the shape, however the 
computer time involved would dramatically increase. The mesh for one period of the structure was 
chosen to be a uniform 85 x 40 x 5 (Figure 8). 



Figure 8 - MAFIA computer mesh simulation of the tunnel ladder structure. 

Computer runs to generate the mesh, set up the eigenvalue matrix, and solve for the 
eigenvalues (frequency resonances) for the two cases, with and without the dielectric rod, were 
completed. In each case for each run the phase shift per period was adjusted from (1/17) -180° to 
(6/17) • 180° to coincide with the experimentally measured frequency resonances. 


COMPARISON OF RESULTS 

Frequency resonance values generated by the MAFIA code were close to experimentally 
measured results. The simulated tunnel ladder structure was then modified to reflect the dimensional 
inconsistencies that were present in the actual structure. The code was run again and frequency values 
were even closer to the experimental results (Figure 9). 
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Figure 9 - Experimental and computer-generated tunnel ladder dispersion characteristics. 

The frequency shifts for both the experimentally measured resonances and the computer- 
generated resonances also were very close and increased with frequency. This is indicative of a 
structure with good interaction impedance for a device where a nominal electron beam voltage may be 
used. The resonance frequencies and shifts in frequencies found using the MAFIA code were slightly 
greater in value than the experimental results. Any discrepancies in results can be attributed to the 
deviation of the actual structure dimensions from the computer model, the computer mesh resolution, 
the dielectric constants for boron nitride and alumina, and the variation in the dielectric cross sections 
along the length of the structure. A much finer resolution would give a better description of the actual 
structure but this is really not needed considering the closeness of the results. Taking into account the 
actual variation and non-uniformity of the experimental tunnel ladder structure dimensions 
demonstrated the feasibility of the MAFIA code for use in the design of these structures. 

CONCLUDING REMARKS 

The ability to perform successful measurements on the actual tunnel ladder structure at the 
operating frequency is a major contribution towards the characterization and analysis of these 
particular slow-wave structures. Previously, scaled versions of the structure were constructed for cold 
testing to determine their suitability to operate successfully at higher frequencies. Also with the 
confidence gained with respect to the computer-generated dispersion characteristics, the ability to 
predict structure characteristics without the need to first fabricate the structure was demonstrated. It 
is much easier, faster, more versatile, and eventually cost effective to model the structure using 
computer-aided design software codes than to fabricate, test, and evaluate a prospective structure 
experimentally to initially determine its potential. Furthermore, variations in the geometry of the 
t structure to achieve the desired characteristics can be incorporated easily and therefore highly efficient 

structures can be designed without prohibitive costs of trial fabrication and testing. 

• It remains for future work to expand the theoretical model of the tunnel ladder structure or 

develop other models to predict with accuracy the dispersion characteristics. Improved fabrication 
methods to reduce structure losses are definitely needed. Although the exact theoretical description of 
many novel structures is extremely difficult, preliminary investigations should be conducted to at least 
approximately characterize the structure in question. 
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With the recent availability of reliable electromagnetic computer programs capable of solving 
Maxwell’s equations in three dimension, a wide variety of periodic slow-wave structures can be designed 
and evaluated easily and without exorbitant cost. The various types of structures that can be 
investigated are limited solely to the imagination of the designer. 
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APPENDIX A 


Mathcad program to find the Mathieu function eigenvalues h and 6 for the elliptical tunnel region. 


Read frequency data file 
frequencies 
number of points N := length (f) 


f := M <0> -10 9 


axial wavenumbers i := 1 .. N 


M ;= READPRN (tlexp) 
co := 2*7C-f 
N = 6 
n i 


p. := — — 

L 17 


speed of light c = 2.99792- 10 8 
permittivity e := 8.854- 10 12 
period L := .0126- .0254 

17 periods of the structure 


Define elliptical parameters. 

semimajor axis semimajor := .03 86-. 0 127 semiminor axis semiminor := .0252- .01 27 

Calculate the radial elliptical coordinate for the ellipse. 


1 , / semimajor + semiminor \ 
pi := —In 1 


2 \ semimaj or - semiminor/ 

Calculate the eccentricity of the ellipse. 


pi =0.78025 


ecc := - 


2 . . 2 
semimaj or - semiminor 


semimajor 

Calculate the interfocal distance. 

d := semimajor- ecc 


ecc =0.75749 


d =3.71335-10“ 


Find the eigenvalues h and b using continued fractions for the TM c01 mode. 


limit := 10 
n := limit.. 1 


m := 0 
s = se 


se := 2 m 
h := 1.9818 


so := 2-m +• 1 
b = 1.52368 


P ratio limit+ 1 := 


16- limit +■ — s 

2 


pratio n 


16-f- s + nj + 2-h 2 - 4-b + h 2 -pratio n+ x 


nratio, 


limit + 1 


16- — s - limit 

\2 


nratio n 


16-(— s-n) + 2-h 2 - 4-b + h 2 -nratk) 


n+1 


Given 


s 2 =b ^2 + pratiOj + nratio^ b := fmd(b) b - 1 


.52368 


Calculate the B 2n coefficients. 


n := 0.. limit 


:= 1 


a n+ i =pra tio n+ i a n B 2 . n := 
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Even solutions of period k 

s - even integer = 2m, allowed values of b = be^ 

See(<(>) = ^B 2 . n cos(2n<t>) 
n 

Jee(h.p) := j^f(- 1)°" m B 2 , n Jn(2n,hcosh(p))| 
n 

Nee(h.p) := J ^][(-l) n - m B 2 . n Yn(2 n > h cosh(p» 
n 

Even solutions of period 2n 

s = odd integer = 2m + 1, allowed values of b = be^ + 1 

Seo(<|») := ^B 2 . n+ |-cos((2 n + 1)<|>) 
n 


Jeo(h.n) 


~J](- 1)" - m B 2n+ ,Jn(2n + 1 ,hcosh(n)) 
n 


Neo(h,|i) := J-^](-l) n ' m B 2 . n+1 Yn(2 iM- l,h cosh(p)) 
n 

Odd solutions of period % 

s = even integer = 2m, allowed values of b = be^ 


Soe(<|>) := y^B 2 . n sin(2 n-(|)) 
n 

Joe(h,|i) := - tanh(p) V|(- l) n_ m (2 n) B 2 n Jn(2 n,h cosh(|i)) 
n 

Noe(h,p) - J^-t anh(p)-^(- l) n ” m (2n)B 2 Yn(2n,hcosh(p)) 
n 

Odd solutions of period 2 tc 

s = odd integer = 2m + 1, allowed values of b = be^ + , 

Soo(<t>) := ]^B 2 . n+ 1 sin((2 n + 1)<|>) 
n 

Joo(h,^i) = - tanh(^) l) n m (2n + l)B 2n+ j-Jn(2*n + l,hcosh(») 

4 2 1 

n 

Noo(h,|i) = -tanh(^)VV l) n " m (2 n + l)B 2 . n + j-Yn^n + 1 ,hcosh(») 

A 2 . * * 
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Condition for E z to be zero on the ellipse. 

Given Jee(h,pl )«0 h = find(h) 


h = 1.9818 


b = 1.52368 


Kp 


,<t>) = d 4 


2 2 
cosh(p) - cos(<t>) 


kl := — 
d 


kc- := 


Field equations for the TM c01 mode. 

E z (|i, 4 >) := Jee(h,p)See(<|>) 

E u (i,p,<t.) = Pj— Jee(h,p)See(<t>) 
dp 


i- 1 1 


H z =0 


H„(i,P,<l>) = ®>j _ j-e-Jee(h,p)— See( 4 >) 

d 4 > 


E^i.P,*)*) : = P;- Jee(h,p) — See(<J>) 

dtp 


Hj/i.n.tlO = fflj_ j-e — Jee(h,p)See(<|>) 
dp 


Power flow through the elliptical cross section. 


Pj = - — — 
1 2 /•- \4 


(S' 


r^i 


0 


r2-7c 


0 


(E (i (i,P.<l ) )'H ( j ) (i,p,<tO + E ^(i.p.^ H^i.p,!)))) dtp dp 


Interaction impedance 
dielectric rod radius 

effective rod semimajor axis a r = 

a r = 3 . 9295 * 10 -4 
b r = 2 . 56537 * 10 -4 


r = . 0125- .0254 
r 


effective rod semiminor axis 


\ r := r -]jl- 


ecc 




ecc 


rod elliptic coordinate p r := - In 

p r = 0.78025 


a r + b r 


l r r/ 


interaction impedance Zj 


C2k 

0 


E (p.tp) 2 l(p,4>) 2 dtp dp 




Plot of interaction impedance based on power flow through the elliptic tunnel region only. 



i 

fj_ ! 1 ° 

Z i 

1 


14 


29.76528 

2 


18.1875 


66.12221 

3 


22.21875 


88.01017 

4 


24.65625 


110.86925 

5 


26.375 


133.2358 

6 


28.09375 


152.53418 


13 



APPENDIX B 


BASIC program to solve for the Mathieu function eigenvalue 6 given eigenvalue h and order s. 


10 DIM PRATIO(15),NRATIO(15),DPRATIO(15),DNRATIO(15),APLUS(15),AMINUS(15), 
BPLUS(15),BMINUS(15) 

20 CLS:INPUT "What is the order s”;S:INPUT "What is h”;H:INPUT "What is the guess for b”;B 
30 INPUT ”What is the number coefficients required”;N 

40 REM-Figure a(n)/a(n-l) 

50 PRATIO(N+l)=-H*H/16/(S/2+N)/(S/2+N) 

60 FOR I=N TO 1 STEP -1 

70 PRATIO(I)=-H*H/(16*(S/2+I)*(S/2+I)+2*H*H-4*B+H*H*PRATIO(I+l)) 

80 NEXT I 

90 APLUS(0)=1:POSASUM=APLUS(0):POSBSUM=0 
100 FOR 1=1 TO N 

110 APLUS(I)=PRATIO(I)*APLUS(I-l) 

120 POSASUM=POSASUM+APLUS(I) 

130 NEXT I 
140 FOR 1=0 TO N 

150 BPLUS(I)=APLUS(I)/POSASUM 
160 POSBSUM=POSBSUM+BPLUS(I) 

170 NEXT I 

180 REM-Figure a(-n)/a(-n+l) 

190 NRATIO(N+l)=-H*H/16/(S/2-N)/(S/2-N) 

200 FOR I=N TO 1 STEP -1 

210 NRATIO(I)=-H*H/(16*(S/2-I)*(S/2-I)+2*H*H-4*B+H*H*NRATIO(I+l)) 

220 NEXT I 

230 AMINUS(0)=1:NEGASUM=AMINUS(0):NEGBSUM=0 
240 FOR 1=1 TO N 

250 AMINUS(I)=NRATIO(I)*AMINUS(I-l) 

260 NEGASUM=NEGASUM+AMINUS(I) 

270 NEXT I 
280 FOR 1=0 TO N 

290 BMINUS(I)=AMINUS(I)/NEGASUM 
300 NEGBSUM=NEGBSUM+BMINUS(I) 

310 NEXT I 

320 REM-Finding b- 

330 DPRATIO(N+1)=0 
340 FOR I=N TO 1 STEP -1 
350 DPRATIO(I)=PRATIO(I)* 

(4-H*H*DPRATIO(I+l))/(16*(S/2+I)*(S/2-|-I)+2*H*H-4*B+H*H*DPRATIO(I+l)) 

360 NEXT I 

370 F=B-H*H/4*(2+PRATIO(l)+NRATIO(l))-S*S 
380 DF=l-H*H/4*(DPRATIO(l)+DNRATIO(l)) 

390 BNEW=B-F/DF:IF ABS(BNEW-B)>.00001 THEN B=BNEW/2+B/2: 

PRINT ”b=”;B:GOTO 40 
400 PRINT ”N”,”a+=”,”B+=”,”a-=”,”B-=” 

410 FOR 1=0 TO N 

420 PRINT I,APLUS(I),BPLUS(I),AMINUS(I),BMINUS(I) 

430 NEXT I 

440 PRINT ”h=”;H,”b=”;B 

450 PRINT "positive a sum=”;POSASUM, "negative a sum=”;NEGASUM 
460 PRINT "positive b sum=”;POSBSUM, "negative b sum=”;NEGBSUM 
470 END 
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APPENDIX C 


Mathieu Eigenfunction Solutions: Periodic in suitable for d real. 

Even solutions about d = 0; b = 6e 2m or be 2m + v 

oo 

Se 2 m (h, cos t?) = 53 B 2n cos(2nt?); E B i n = 1 

~ n 

n = 0 

oo 

5e 2ro + X (A, cos t?) = 53 B 2 n + 1 cos[(2n + l)t?]; E 5 2 „ + 1 = 1 

n = 0 " 

Odd solutions about d = 0; b = &o 2m or ^°2m + 1* 

oo 

So 2m (h , cos tf) = 53 B 2n sin(2m?); E( 2 »)*2n = 1 

n = 1 
oo 

S°2 m + l(^» cos = ^ v ^2 n + 1 s ^ n [(2^ + l)^]i (^ n “b ^)^2n + 1 “ * 

n = 0 n 

Corresponding Radial Solutions 

For /i = id and for values of the B ' s and of 6 corresponding to the angular functions Se , 5o. 

Even functions: 

Je 2m (h, cosh /i) = fi/2 £ ( - cosh /i) 

n = 0 
oo 

e 2m + l(^> C0S ^ /*) = '\' k /2 ( ” 1)” m ^2n + 1^2n + l(^ C0S ^ /*) 

n = 0 

Odd functions: 

Jo 2m {h, cosh /i) = ^7r/2 tanh ^53 ( “ 1)" “ m (2n)2? 2B J 2n (fc cosh A*) 

n = 1 


•A °2m + l( A ’ cosh A 1 ) = tanh A* 53 ( " 1 )" m ( 2n + + l J 2n + l( h cosh A*) 

n = 0 
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